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ABSTRACT 

We propose that the rotational kinematics of the globular cluster system (GCS) 
in M31 can result from a past major merger event that could have formed its bulge 
component. We numerically investigate kinematical properties of globular clusters 
(GCs) in remnants of galaxy mergers between two disks with GCs in both their disk 
and halo components. We find that the GCS formed during major merging can show 
strongly rotational kinematics with the maximum rotational velocities of ~ 140 — 170 
km s~^ for a certain range of orbital parameters of merging. We also find that a rotating 
stellar bar, which can be morphologically identified as a boxy bulge if seen edge-on, 
can be formed in models for which the GCSs show strongly rotational kinematics. We 
thus suggest that the observed rotational kinematics of GCs with different metallicities 
in M31 can be closely associated with the ancient major merger event. We discuss 
whether the formation of the rotating bulge/bar in M31 can be due to the ancient 
merger. 

Key words: Magellanic Clouds - galaxiesistructure - galaxies:kinematics and dy- 
namics - galaxies :halos - galaxies:star clusters 



1 INTRODUCTION 

The structural and kinematical properties of GCSs are con- 
sidered to provide valuable information on the formation and 
evolution of their host galaxies (e.g., Searle & Zinn 1978; Ro- 
manowsky et al. 2009) . Recent observational studies on kine- 
matical properties of GCSs in galaxies and their comparison 
with theoretical and numerical works have advanced our un- 
derstanding on mass distributions of galaxies (e.g., Pierce et 
al. 2006; Romanowsky et al. 2009), formation of elliptical 
galaxies (e.g., Bekki et al. 2005), and formation and evolu- 
tion of dwarf galaxies (e.g., Beasley et al. 2009). A growing 
number of observational data sets on kinematical properties 
of GCSs have been now accumulated not only for galaxies 
in the local group but also for nearby galaxies (Brodie & 
Strader 2006 for a recent review) so that we can discuss 
formation processes of their host galaxies and their depen- 
dences on galactic global properties (e.g., Hubble types) in 
more detail based on the observations. 

One of the intriguing properties of GCSs in galaxies 
is the observed rotational kinematics of the GCS in M31 
(e.g., Huchra et al. 1982; Perrett et al. 2002). Recent 
wide-field surveys of the GCS in M31 (Lee et al. 2008) 
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have confirmed that the GCS composed both of metal-poor 
([Fe/H] < —0.9 in their criterion) and metal-rich (Fe/H] 
> —0.9) objects has a rotational amplitude of ~ 190 km 
s^^. Such kinematics are in a striking contrast with the low 
rotation of the Galaxy's GCS (e.g., Armandroff 1989) with 
an origin that remains unclear. 



The purpose of this paper is to discuss (i) why the GCS 
of M31 shows such a large amount of rotation and (ii) what 
implications it has for the formation and evolution of M31. 
Previous numerical works showed that rotational kinemat- 
ics of GCSs in elliptical galaxies can be due to past major 
merger events that formed the galaxies (e.g., Hernquist & 
Bolte 1993; Bekki et al. 2005). We thus consider that the 
rotational kinematics of the GCS is closely associated with 
a major merger event in M31 long ago and thereby inves- 
tigate whether galaxy merging can reproduce well the ob- 
served kinematics of the GCS. We also discuss whether the 
observed rotating bulge/bar in the M31 (e.g., Beaton et al. 
2007) can be due to the major merger event responsible for 
the formation of the GCS with rotational kinematics. 
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Table 1. The values of model parameters and a brief summary of 
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136 123 the standard model 

114 133 

95 92 

107 38 unequal-mass merger 
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The mass ratio of dark matter halo to stellar disk in a galaxy. 
^ The mass ratio of two disks in a galaxy merger. 

The pericenter distance of a merger in units of the disk size . 

The orbital eccentricity of a merger. 
^ "PP" and "PR" represent prograde-prograde and prograde-retrograde merging, respectively. 
f The maximum rotational velocity of DGCs (GCs initially in disks) in a merger remnant. 
9 The maximum rotational velocity of HGCs (GCs initially in halos) in a merger remnant. 



2 THE MODEL 

Since the numerical methods and techniques we employ for 
modeling dynamical evolution of mergers between two disks 
with GCs have already been detailed elsewhere (Bekki et al. 
2005; Bekki & Forbes 2006), we give only a brief review here. 
The progenitor disk galaxies that take part in a merger are 
given a dark halo, a thin exponential disk, and GCs initially 
in disks (referred to as "DGCs" from now on) and in halos 
("HGCs"). The total mass and size of an exponential disk 
with no bulge are Md and Rd, respectively. 

We consider that the total stellar mass of a merger rem- 
nant should be similar to the total mass (Mb) of the present 
bulge of M31 (Mb = 3 - 4 x lO^^M©; e.g., Seigar et al. 2008; 
Geehan et al. 2006). We thus adopt Md = 2 x 1O^°M0 as 
a reasonable value in the present study. We adopt the den- 
sity distribution of the NFW halo (Navarro, Frenk & White 
1996) with the mass of Mdm and a concentration parame- 
ter determined by Mdm according to the formula derived by 
recent CDM simulations (e.g., Neto et al 2007). The stellar 
disk has the scale length (Ro) of 0.2i?d. 

The present bulge-less disk model would be reasonable 
for less massive disk galaxies with no/little bulges like the 
Large Magellanic Cloud and Magellanic-type galaxies with 
total masses of ~ IO^^Mq, but it would not be so realistic 
for luminous disk galaxies like the merger progenitor disks of 
M31. We however conjecture since the bulge masses are quite 
small in comparison with total masses of galaxies (inclusive 
of dark matter halos), final kinematics of GCSs in merger 
remnants would not depend strongly on whether or not we 
include bulges in initial disks. Thus, as long as we mainly 
discuss global kinematical properties of GCSs in merger rem- 
nants, the present model can be regarded as an appropriate 
one. 

Previous observations revealed that less luminous disk 
galaxies, which are considered to form the bulge component 
of M31 with rotational kinematics of the GCS in the present 
study, have disky distributions with rotational kinematics in 
their GCs (e.g.. Freeman et al. 1983; Olsen et al. 2004). We 
thus consider the presence of DGCs in merger progenitor 
disks in the present study: DGCs were not considered in 
our previous works (e.g., Bekki et al. 2003a; 2005) and such 
DGCs with strongly rotational kinematics are not observed 
in the metal-rich GCs of the Galaxy. The DGCs have the 
same exponential distribution and rotational kinematics as 



field stars in the stellar disk of their host galaxy. The initial 
rotational amplitude of DGCs is ~ 170 km s~^ for Md = 
2 X IO^^Mq and Mdm/A/d = 9. 

The Galactic HGCs and the stellar halo have similar 
radial density profiles of p(r) cc r~^'^ (van den Bergh 2000). 
We therefore assume that the HGCs in our galaxies have 
a power-law profile with an exponent of —3.5 and a half- 
number radius of 1.47?o (which is ~ 5 kpc for the Galaxy 
and thus consistent with observations). The HGCs are as- 
sumed to have isotropic velocity dispersions determined by 
the mass distribution of the galaxy. Total numbers of DGCs 
and HGCs in a galaxy are set to be 100 and 100, respectively. 

The mass ratio of the two disks (m2) in a merger is as- 
sumed to be a free parameter. In all of the simulations of pair 
mergers, the orbit of the two disks is set to be initially in the 
xy plane and the distance between the center of mass of the 
two disks is assumed to be 12Rd. The pericenter distance 
(vp) and the eccentricity (cp) in a pair merger are assumed 
to be free parameters that control orbital energy and angu- 
lar momentum of the merger. The spin of each galaxy in a 
merger is specified by two angles 9i and ^i, where suffix i 
is used to identify each galaxy. 8i is the angle between the 
z axis and the vector of the angular momentum of a disk. 
(fji is the azimuthal angle measured from the x axis to the 
projection of the angular momentum vector of a disk onto 
the xy plane. 

We mainly show the results of "the standard model" 
with Md = 2 X 1O^°M0, Mdm/Md = 9, m2 = 1, rp = 2i?d, 
Bp = 0.72, = 30°, 6*2 = 45°, (pi = 45°, and 02 = 120°. 
We describe the results of this standard model with an or- 
bital configuration similar to "prograde-prograde" merging 
(in which the orbital spin axis of the merger is parallel to 
intrinsic spin axes of the two disks), mainly because the fi- 
nal GCS can clearly show strongly rotational kinematics as 
observed in M31. We also show the results of the model 
with a "prograde-retrograde" orbital configuration in which 
02 = 225° and other parameters are exactly the same as 
those in the standard model. 

In order to discuss the origin of the central bar/bulge of 
M31, we investigate in which models the merger remnants 
show rotating stellar bars. Although we investigate 34 mod- 
els, we show the results of five representative models which 
either show clearly rotating bars or have no such rotating 
components in merger remnants: comparison between these 
models enables us to understand better the role of an an- 
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Figure 1. Distributions of disk GCs (DGC, red) and halo ones 
(HGC, blue) of a galaxy merger projected onto the x-y plane (i.e., 
orbital plane of the merger) for larger (upper) and smaller (lower) 
scales of view in the standard model. The time T shown in the 
upper left corner of each panel indicates the time that elapsed 
since the simulation starts. Green circles indicate 20 kpc from 
the mass center of the merger. 



cient merger event in the formation of the bulge component 
in M31. 

The model parameters and some brief results are given 
in the Table 1. The total particle number for a major merger 
in a simulation is 10^ and the simulation is carried out on the 
latest version of GRAPE (GRavity PipE, GRAPE-7) which 
is the special-purpose computer for gravitational dynamics 
(Sugimoto et al. 1990). In estimating the GCS kinematics, 
the merger remnant is viewed near to edge-on, and binned 
major-axis profiles are constructed. In order to have enough 
objects in each bin, GCs at all minor axis distances are in- 
cluded, which also approximately replicates the comparison 




Figure 2. Radial profiles of line-of-sight velocities (V, solid) and 
velocity dispersions (cr, dotted) for DGCs (red) and HGCs (blue) 
for the merger remnant at T = 4.2 Gyr in the standard model. 
The dotted horizontal black line represents V (cr) = km . 
The observed v^ot ~ 140 km s"'^ for the full sample of GCs in 
M31 (Perrett et al. 2002) is shown by a dashed green line for 
comparison: Lee et al (2008) show significantly larger Vrot (~ 188 
km s-l) of the GCS for \Y\ ^ kpc. The locations of the filled 
circles indicate where the observed velocity profile becomes flat 
and the error bars are observational ones (±13 km s~^). 
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Figure 3. The same as Fig. 2 but for for dark matter (red) 
and disk stars (blue). The bulge observational results are from 
McElroy (1983) and the locations of the filled circles indicate the 
observed outermost region {R = 9.76 ) at the position angle of 
45°. 



observations. Formation of metal-rich GCs through dissipa- 
tive gas dynamics of galaxy merging (e.g., Bekki et al. 2002) 
and adiabatic contraction of the GCS in M31 due to later 
massive growth of the disk component could increase appre- 
ciably the rotational amplitude V of the GCS. The quanti- 
tative estimation of this effect is not done here and will be 
done in our future studies. It should be stressed that this 
paper is meant to be schematic rather than reproducing the 
observed galaxy in a fully self-consistent manner. 



3 RESULTS 

Fig. 1 shows that both DGCs and HGCs can be spatially 
mixed to form a new GCS during violent dynamical rela- 
tion that results in transformation from two disks into one 
spheroidal galaxy. Owing to angular momentum redistribu- 
tion during major merging, not only HGCs but also DGCs 
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Figure 4. The same as Fig. 2 but for tlie model witli m2 = 0.5. 
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Figure 5. Distributions of stars projected onto tlie x-y (left) and 
x-z (right) planes at T = 4.2 Gyr in the standard model. Only 
one for every ten particles is shown for clarity. 



can be transferred to the outer halo region and thus seen 
there in the merger remnant. The half- number radii for the 
DGCs and the HGCs are both ~ 5.0 kpc in the merger rem- 
nant, which is consistent with the observations by Battistini 
et al. (1993). As the stellar remnant shows a rotating bar 
(discussed later in §4), DGCs and HGCs also show barred 
structures, appreciably flattened shapes (qualitatively con- 
sistent with observations), and figure rotation. 

Fig. 2 clearly shows global rotation both in DGCs and 
HGCs, though the radial profiles of V and a do not smoothly 
change owing to the small numbers of GCs at each bin. The 
maximum V (Fmax) and a (cTmax) for DGCs (HGCs) within 
the central 16 kpc are 136 km s~^ (123 km s"'^) and 102 
km s~^ (106 km s~^), respectively. Therefore Vmax/o"max for 
DGCs and HGCs are 1.3 and 1.2, respectively, which means 
that the final GCS has strongly rotational kinematics. The 
velocity profile becomes fiat {V ~ 100 km s~^) at R ^ 4 kpc. 
The simulated two-dimensional line-of-sight velocity map of 
the GCS shows clearly global rotation. 

We here compare the simulated rotational kinematics 
with the observed one in M31 (e.g., V — 138 ± 13 km s~^ in 
Perrett et al. 2002) . Observational error bars in V axe not so 
small and V ranges from 98 km s~^ to 188 km s"'^ for Okpc 
^ \Y\ ^ 5kpc (Lee et al. 2008), where \Y\ is the projected 
vertical distance from the M31 disk plane. The simulated 
rational amplitude of the GCS is slightly smaller than the 
observed one by Perrett et al. (2002): It should be noted 
that the observational results depend on the details of how 
the data are binned and fitted. 

We consider that the best model needs to reproduce the 
above V of ~ 140 km s~^ Fig. 3 shows that even the dark 



matter halo of the merger remnant can have a small amount 
of rotation {V ~ 34 km s~^) owing to the redistribution of 
angular momentum (i.e., conversion of orbital angular mo- 
mentum of merging two galaxies into internal one of the 
merger remnant). This result suggests that if the two galax- 
ies have hot gaseous halos, then the remnant is highly likely 
to have a slowly rotating gaseous halo. Fig. 3 also shows that 
the stellar component of the merger remnant has a signifi- 
cant amount of rotation (Vmax ~ 100 km s^^), which means 
that the spheroid is rapidly rotating (i.e., rotating bulge is 
formed from major merging). 

The reason for the smaller V of the dark matter halo 
in comparison with the GCS is due largely to the difference 
in the initial spatial distribution between the halo and the 
GCS, demonstrating that GC rotation would not be used to 
infer dark matter halo rotation. The inner part of the merger 
can more strongly spin-up during and after major merging 
owing to (i) the prograde-prograde orbital configuration of 
the merging and (ii) the development of the rotating bar. 
The distribution of the GCS in the merger progenitor disk 
is by a factor of 4 more compact that that of the halo so 
that the GCS can more strongly spin-up: most of the in- 
dividual GCs can obtain a larger amount of intrinsic spin 
angular momentum with respect to the center of the merger 
remnant. 

Fig. 4 shows that strongly rotational kinematics of GCs 
can be seen only in DGCs in the model with — 0.5 
(Model 4): Vmax is 107 km s^^ for DGCs and 38 km s"^ for 
HGCs. The higher and lower Vmax in DGCs and HGCs, re- 
spectively, are confirmed in other unequal-mass merger mod- 
els (e.g., 7712 = 0.3). It should be noted here that the model 
1 shows only slightly higher Vmax in DGCs, which is due 
to different kinematics between DGCs and HGCs (i.e., im- 
tial global rotation only in DGCs). These results therefore 
mean that if metal-poor and metal-rich GCs originate from 
HGCs and DGCs, respectively, then kinematical properties 
of GCs in the remnants of galaxy merging with smaller m2 
(or unequal-mass merging) can be significantly different be- 
tween metal-poor and metal-rich ones. These results further- 
more imply that the observed rotational kinematics both 
for metal-poor and metal-rich GCs in M31 can give some 
constraints on the mass-ratio (m2) of two disks in galaxy 
merging that could have occurred in M31. 

The model 1 does not explain well the observed to- 
tal halo mass: later numerous accretion events of dwarfs 
with few GCs and stars are required to increase signifi- 
cantly the total mass after the merging. The models with 
larger Mdm/Afd can show higher Vmax in DGCs owing to 
the initially higher circular velocities of the stellar disks. 
For example, the model with Afdm/Md = 19 (Model 5) 
shows Vmax = 174 km s^^ and amax = 134 km s^^: it 
should be stressed here that the total mass of the remnant 
IS 8x 10"Mq thus can be more consistent with observations 
than Model 1 in terms of the total mass of M31. A possi- 
ble reason for the lower Vmax (91 km s~^) for HGCs in the 
model 5 is that more strongly self-gravitating stellar disks 
can also play a role in increasing global rotation of HGCs: 
Such a role is weaker in the model 5 in which the disk is 
much more weakly self-gravitating. 

Models with different orbital configurations can show 
stronger rotational kinematics in DGCs and HGCs, if larger 
Tp are adopted. For example, the model with a prograde- 
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retrograde orbital configuration (Model 2) shows Vmax = 114 
km s~^ and amax = 102 km s~^ but does not show a bar. The 
models with smaller rp (e.g., Model 3) show smaller Knax in 
DGCs, which implies that the observed Vmax can give some 
constraints on the orbital parameters for galaxy merging 
that occurred in M31. Thus only the models with larger r-p 
can better reproduce the observed kinematical properties of 
the M31 GCS. 



4 DISCUSSION AND CONCLUSIONS 

Fig. 5 shows that the stellar remnant looks like a bar if it 
is viewed from face-on in the standard model of the present 
study: the bar is confirmed to have figure rotation. The stel- 
lar distribution viewed from edge-on appears to have a flat- 
tened spherical body, which can be identified as a bulge. The 
present study thus implies that M31's observed bulge/bar 
(e.g., Beaton et al. 2007) can be formed from an ancient 
major merger event. It should be noted here that dissipativc 
gas dynamics which can determine the final morphological 
properties of merger remnants (e.g., boxy or disky shapes) 
is not included in the present study. 

If the inner bar/bulge of M31 was really formed from 
ancient major merging before its disk formation, then the 
later development of the stellar and gaseous disk may well be 
significantly infiuenced by dynamical action of the already 
formed bar. Also, later slow gas accretion and the resultant 
disk formation in M31 could change structural and kinemat- 
ics properties of the already developed GCS to some extent: 
it would be possible that adiabatic compression of the GCS 
by later gradual development of the disk can enhance the 
rotational amplitude of the GCS. It is our future study to 
numerically investigate disk formation and evolution of M31 
under the presence of the already formed bar. 

The present work suggests that the hot diffuse halo 
gas recently detected by Chandra (Li & Wang 2007) can 
have a significant amount of rotation resulting from angu- 
lar momentum redistribution of the possible major merger. 
Furthermore, the observed extensive HI cloud population of 
M31 (e.g., Thilker et al. 2004; Westmeier et al. 2007) can also 
have rotational kinematics if they originate from stripped 
HI gas from the merger progenitor disks. Given that ma- 
jor merging can form very extended stellar halos (Bekki et 
al. 2003b; Bekki & Peng 2006), the observed extended stel- 
lar halo in M31 (e.g., Ibata ct al. 2007) would have fossil 
information in the possible ancient major merger event. 

If M31 was formed from multiple and sequential merg- 
ing of dwarf satellites from random directions, it seems to 
be highly unlikely that the final GCS has a large amount 
of global rotation. Therefore, if the rotational kinematics 
of the GCS derived from previous observations is further 
confirmed by ongoing observational studies, it then suggests 
that a dramatic physical process is responsible for the ob- 
served globally organized motion of GCs. Our simulations 
imply that an ancient major merger of two disks with GCs 
with different metallicities could be responsible for the enig- 
matic kinematics of the GCS in M31, as well as the inner 
bar or bulge. 
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